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ABSTRACT 

Four microlensing collaborations are presently searching for compact matter in the 
Galaxy and all have detected possible candidates. Using the detection efficiencies re- 
cently published by the MACHO and OGLE collaborations, we present Monte-Carlo 
calculations of the expected optical depth, rates and timescales, along the LMC and 
Galactic bulge lines of sight, for dark matter in a four-component 'standard Galaxy' 
model with a spherically-symmetric halo and spheroid. Using the typically observed 
event durations we show that, whilst the halo fraction comprised of compact matter 
is likely to be fh < 0.4, a 'no halo compact matter' hypothesis is ruled out at greater 
than the 80% confidence level, unless the LMC itself has a substantial halo of such 
objects. On the basis of the timescales observed by OGLE towards the bulge we find 
the rate predicted by the model to be in good agreement with the number of OGLE 
detections. We compute lens mass probability distributions for the various components 
and compare these estimates with current observational and theoretical constraints on 
the mass scale of baryonic dark matter. We assess the uniformity of the amplification 
distributions for the published EROS, MACHO and OGLE events and find that they 
are quite consistent with the microlensing hypothesis, although the OGLE candidate 
selection criteria mean that its data are particularly sensitive to photometric selection 
effects. The EROS team has recently placed strong limits on the density contribution 
of very low mass halo objects from their short timescale CCD search. On the basis of 
a recent study of the fiux amplification of finite-size sources by Simmons, Newsam & 
Willis (1995) we suggest that EROS may have detected up to 5 low-amplification events 
due to halo lenses with mass m ~ 10~^ Mq. 
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INTRODUCTION 



The dark matter problem in our own and other galaxies is 
now well established; e.g. Rubin & Ford (1970), Ostriker & 
Peebles (1974), Ashman (1992). The inferred mass of galac- 
tic dark matter is typically estimated to be an order of mag- 
nitude more than that of the luminous parts of galaxies. 
Many possible solutions have been invoked from the mini- 
mal assumption of dark baryonic matter (Carr 1994) to the 
existence of a new form of weakly interacting matter [often 
referred to as cold dark matter (CDM) or weakly interacting 
massive particles (WIMPs)]. 

It is known from the observed abundances of the light el- 
ements that baryonic matter cannot provide more than 10% 
of the closure density of the Universe (Walker et al. 1991), 
but this still admits the possibility that it could provide all 
or some of the galactic dark matter. Indeed, the same obser- 
vations imply a lower limit for the baryon density which is 
three times greater than the density in visible baryons (Per- 
sic & Salucci 1992) and so require some of the dark matter 
to be in baryonic form. Tentative measurements of a high 



deuterium abundance in a high-redshift primordial gas cloud 
by Songaila et al. (1994) and Carswell et al. (1994) imply a 
baryon density which is comparable to the visible density, 
but the discrepancy with solar system abundances implies 
that these observations may be due to an intervening hydro- 
gen cloud at lower redshift (Steigman 1994). Subsequently, 
observations of another high-redshift cloud have inferred a 
much lower deuterium abundance (Carswell 1994). A com- 
parison of the amount of gas present in some galaxy clusters 
with predictions from numerical simulations implies that the 
contribution of baryons to the total cluster mass may exceed 
the amount compatible with standard nucleosynthesis calcu- 
lations (White et al. 1993), unless the density of the universe 
is only around 20% of the closure density. It therefore seems 
likely that baryonic matter is relevant on galactic scales. 

The nature of such baryonic matter is unknown, but 
there are a number of observational and theoretical con- 
straints on many of the possible candidates; e.g. Carr (1994), 
Dalcanton et al. (1994), Silk (1994). Paczyhski (1986) sug- 
gested that microlensing could be used to look for compact 
dark matter in the Galactic halo by simultaneously mon- 
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itoring the brightness of 10^"'^ stars for a few years for 
signs of variations indicative of the microlensing effect. Four 
experiments are now in progress: MACHO (Alcock et al. 
f993); EROS (Aubourg et al. f993); OGLE (Udalski et al. 
f993); and DUO (D. Valls-Gabaud, private communication). 
To date these teams have detected more than 70 candidate 
events. 

The plan of this paper is as follows. In Section 2 we 
discuss the lensing properties of dark matter in a 'stan- 
dard Galaxy' model comprising a spherically-symmetric 
halo, thin and thick discs and a spherically-symmetric 
spheroid/bulge component. We use the recent efficiency esti- 
mates from the MACHO and OGLE collaborations and the 
typically observed event durations to compare the predic- 
tions of this model to the observations, and assess to what 
extent the observations disagree with the model predictions. 
We also discuss the implications of the current results from 
EROS short-timescale CCD search which is sensitive to very 
low mass objects. In Section 3 we calculate mass probability 
distributions for the lenses along both the LMC and bulge 
lines of sight and discuss their consistency with other theo- 
retical and observational constraints on the mass scale bary- 
onic dark matter. We also test the distribution of event am- 
plifications observed by EROS, MACHO and OGLE against 
the expected distribution for microlensing events. Section 4 
discusses the results and implications of the previous sec- 
tions. 



2 OPTICAL DEPTH AND LENSING RATE 
TOWARDS THE LMC AND BULGE 

The number of lensing events in progress at a particular 
instant in time is given by the optical depth: 



TT f Rl(m,x,L) , , (i^TV, 



dxdldbdL, (2.1) 



where the integral sign implies integration over the source 
distance L (measured from the observer), lens distance 
X < L, Galactic longitude I and latitude b. pi denotes the 
lens density, m the lens mass, Ns the number of sources mon- 
itored and = [4Gmx(i - x)/c^i]^/^ the Einst ein radius, 
which is the radius of the image ring formed (as measured 
along the lens plane) when the lens is perfectly aligned along 
the observer-source line of sight. If the lens and source are 
separated by a distance uRe then the source brightness is 
amplified by a factor 



A 



m(m2 -^4)1/2 



(2.2) 



Re provides a natural lensing length scale and from here on 
we shall define an event to be where the lens-source separa- 
tion is less than Re [i.e. m < f which, from eqn (2.2), implies 
A > f.34]. 

The number of events due to lenses of mass m per unit 
time interval is the lensing rate: 

r = Tr / ' p,(x,l,b)(Vt(x,L))e 



dl db dL 



dx dl db dL, 



(2.3) 



where Vt is the relative tangential velocity of the lens with 
respect to the line of sight. Vt depends on the respective 
motions of the lens and source and their relative orienta- 
tion 9 (Griest f99f ; Kiraga & Paczyhski f994) and so there 
is an implied summation over all possible source and lens 
velocities, as well as over I, b, L and x < L. The average 
event timescale is (te) = (Re/Vt) = 2t / irT and is thus de- 
termined by the optical depth and lensing rate. Note that 
we are here defining the event timescale as the Einstein- 
radius crossing time. The MACHO collaboration adopt the 
Einstein-diameter crossing time as their definition and so all 
timescales in this paper need to be multiplied by 2 before 
being compared to the MACHO results. 

There are many possible locations for the lensing events. 
Towards the LMC possible sources include the Galactic halo 
(Paczyhski f986; Griest f99f; De Riijula, Jetzer & Masso 
f99f); the thin or thick discs (Paczyhski f99f; Griest et al. 
f99f ; Gould f994a; De Riijula, Giudice, MoUerach & Roulet 
f994); the spheroid (Giudice, MoUerach & Roulet f994); and 
the LMC itself (Gould f 993; Sahu f 994; De Riijula, Giudice, 
MoUerach & Roulet f994). The halo, disc and spheroid may 
also contribute to lensing towards the Galactic bulge, along 
with the bulge itself (Kiraga & Paczyhski f 994; Evans f 994; 
De Riijula, Giudice, MoUerach & Roulet f994). 

Gould, Miralda-Escude & Bahcall (f994) have investi- 
gated the relative contributions of the halo, thin disc and 
thick disc components to the optical depth towards the 
LMC, SMC and bulge. A similar analysis of the lensing rates 
could in principle provide a more critical comparison, for a 
particular Galactic model, since, not only does the rate de- 
pend on the lens density distribution, it also depends on the 
lens mass distribution and the lens and source velocity pro- 
files. However, a naive analysis which assumes that events 
of all timescales are uniformly sampled by the microlensing 
experiments is likely to be misleading since, in practice, sam- 
pling and photometric selection effects mean that microlens- 
ing observations do not uniformly sample the lensing rate. 
The frequency of individual observations dictates the mini- 
mum observable timescale te(min), whilst the total duration 
of the experiment determines the maximum observable event 
duration te(max). Timescales between te(min) and te(max) 
are sampled with an efficiency e(te) < 1 (Griest f99f). Pho- 
tometric effects, such as poor seeing or excessive 'crowding' 
of the stellar images, also mean that some event durations 
will be observed in preference to others. 

The different density and velocity profiles of the various 
components give rise to different characteristic lensing event 
durations, and so the effect of the sampling and photomet- 
ric efficiencies can be to enhance the observed rate of lenses 
from one component relative to another. The MACHO and 
OGLE collaborations have recently published their detection 
efficiencies (Bennett et al. f995; Alcock et al. f995; Udalski 
et al. f994) and so we shall use these in examining the rela- 
tive contributions of a four-component Galactic model to the 
observed rate and optical depth along the LMC and bulge 
lines of sight. 

In evaluating the lensing rate towards the LMC we will 
neglect the contribution due to lenses in the LMC itself and 
consider only Galactic sources, although we shall discuss the 
likely LMC contribution. In the case of lensing towards the 
bulge we follow Kiraga & Paczyhski (f 994) in taking account 
of the line-of-sight depth of the bulge sources. 



Galactic microlensing 3 



2.1 Galactic model 

The three main components comprising our Galaxy are 
the rotationally-supported disc, the partially rotationally- 
supported visible halo and the dark outer halo. There is 
some evidence that the disc comprises both a 'thin' and 
'thick' component, whilst the visible halo can be sub-divided 
into a central bulge surrounded by a spheroidal component. 
[For reviews see Bahcall (1986) and Gilmore (1989).] We 
shall utilise simple analytic models for these components. 
For the dark halo we shall adopt the softened-isothermal 
sphere of Bahcall & Soneira (1980): 



/9(halo) = po 



a\ + Rl 



(2.4) 



where po = 0.01 Mq pc"^ is the local density, Ro = 8 kpc 
is the Sun's Galactocentric distance and au — 5 kpc is the 
halo core radius. We assume that the halo lenses have an 
isotropic random velocity distribution (as measured in the 
rest frame of the Galaxy) which, in two dimensions, resem- 
bles the Rayleigh distribution: 



P(V,)dV, = 2 (V,/Vc) exp[-(V,/Vcf] (dV,/Vc), 



(2.5) 



where Vc = 220 km s~^. The optical depth and lensing rate 
produced by a halo of the form of eqns (2.4) and (2.5) has 
been calculated by Griest (1991) and has recently been con- 
strained by the MACHO collaboration (Alcock et al. 1995). 
Here, we shall take account of other possible sources of lens- 
ing in assessing the contribution of the halo. 

For the thin and thick discs we take the density distri- 
bution to be 



p(disc) = Po 



X exp 



Ro — (x^ cos^ b — 2xRo cos fccos I -\- R^Y^^ 



X sin I b I 



H 



(2.6) 



where h and H are the disc scale length and scale height, 
respectively, and po = Pvis + Pdark, with p,j,s the contribu- 
tion of the visible matter and pdark the contribution of the 
dark matter. For the thin disc we take po = 0.2 Mq pc~'^, 
with Pdark = 0.1 Mq pc~'^ and H = 0.3 kpc, whilst for 
the thick disc we assume po = 4 x 10"'^ Mq pc~'^, with 
pdarfc = 2 X 10~^ Mq pc~^ and H = 1.0 kpc. We adopt 
h = 3.5 kpc for both components. We assume that both com- 
ponents are rotationally supported with Vrot = 220 km s~^ 
for the thin disc and Vrot = 170 km s~^ for the thick compo- 
nent. In addition, we take both discs to have a random veloc- 
ity component of the form of eqn (2.5) with Vc = 30 km s~^ 
and 60 km s~^ for the thin and thick disc, respectively. 

For the spheroid and bulge (which, form here onwards, 
will be collectively referred to as the bulge) we adopt the 
simple spherically-symmetric model used by De Riijula et 
al. (1994) where 



p(bulg 



Po 



1/2 



+ ^0 



1/2 



1/2 I / 2 
flfc + [X^ 



2xRo cos 6 cos I -I- Riy/^ 



(2.7) 



with Po = Pvis + Pdark- There is some uncertainty as 
to the amount of dark matter in the bulge and so we 



shall discuss two scenarios: po = 10~ Mq pc~ , with 
Pdark = 8.75 X 10~* Mq pc~'^ (which, from here on, will 
be referred to as the 'dark matter-dominated' bulge), and 
Po = 2.5x10"* Mq pc"^ with pdark = 1.25x10"* Mq pc"^ 
The core radius is taken to be at = 0.17 kpc. We take the 
bulge to have a rotational motion Vrot = 100 km s"^ and a 
random component given by eqn (2.5) with Vc = 100 km s"^. 
Spherical symmetry is likely to be an over simplistic as- 
sumption since much of the observational evidence points 
towards a flattened or perhaps bar-like bulge [e.g. Blitz & 
Spergel (1991); Weiland et al. (1994)]. However, calculations 
based on spherical-symmetry at least provide a 'bench mark' 
against which any asymmetry can be quantified. One of the 
aims of this paper is to evaluate to what extent microlensing 
observations provide evidence for bulge asymmetry. 

The visible component of the local spheroid density is 
normalised to l/800th that of the thin disc, and the visible 
component of the thick disc density to l/50th. The combined 
column density (dark and visible matter) of the four com- 
ponents within 1.1 kpc of the Galactic plane is 83 Mq pc"^, 
consistent with observational constraints (Bahcall, Flynn & 
Gould 1992). The total column density in the thin and thick 
discs is 64 Mq pc"^, which is also consistent with rotation 
curve constraints. We shall refer to this as the 'standard 
model', since it is very much along the lines of the Galac- 
tic models advocated by Bahcall & Soneira (1980) and by 
Kuijken & Gilmore (1989). The main uncertainty regarding 
the standard model concerns the existence of the thick disc. 
However, it turns out that the lensing contribution of the 
thick disc is anyway rather small. 

We are concerned here with the lensing properties of 
the dark matter components of this model. The visible com- 
ponents also make some contribution to the lensing rate, 
although the number of achromatic lensing events (which 
generally require the lens to be invisible) is small. This fact is 
reinforced by recent limits on the lensing contribution of M- 
dwarfs (Bahcall, Flynn, Gould & Kirhakos 1994; Hu, Huang, 
Gilmore & Cowie 1994). 



2.2 Source distribution 

As well as specifying the lens distribution, we also need 
to do the same for the sources. In the case of lensing to- 
wards the LMC, the sources can be assumed to be all at the 
same distance, since the distance to the LMC is significantly 
greater than its line-of-sight depth. The angular distribution 
of the sources is also unimportant. We shall therefore adopt 
i = 50 kpc, I = 280° and b = -33° for the LMC sources. 
We shall, however, follow Griest (1991) in taking account of 
the tangential motion of the LMC. 

In calculating the optical depth and lensing rate towards 
the bulge, the line-of-sight depth can provide an important 
correction (Kiraga & Paczyhski 1994). We assume that all 
sources reside in the bulge although one should also expect 
a small contribution to the lensing rate from sources in the 
disc. Taking the number density of sources per unit V-band 
luminosity to be dng/dC oc between Ci and £„, the 

variation of the number of sources with I, b and L is given 
by 

d^Ns oc [^i"" - C];;f„(L)\ps(l,b,L)L^ cosbdldbdL 

(/?/l), (2.8) 
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where the source density ps is given by eqn (2.7) and 
i^mtn = max[£i, A'(i/kpc)^] is the minimum observable lu- 
minosity, with K = io5 9-o *^('"")/:q (neglecting extinc- 
tion) and V{Iim) ~ +21 mag the limiting apparent V-band 
magnitude of the observation. For our calculations we adopt 
/? = 2. Since Crmn ^ for all reasonable values of L and 
Cu, eqn (2.8) turns out to be insensitive to for /? > 1. 
Furthermore, assuming Ci < 10"'^ Cq (corresponding to a 
source mass nis < 0.1 Mq), Crmn > Ci everywhere along the 
line of sight except very close to the observer (L < 0.5 kpc). 
This means that eqn (2.8) is also insensitive to the precise 
value of Ci and so d^Ns oc ps(L)L^'^^-'^^ = ps(L). 

2.3 Results 

To evaluate eqns (2.1) and (2.3) we used Monte-Carlo in- 
tegration. The (Gaussian) error estimate was typically less 
than 1% for 10^ evaluations. In reality, the error distribu- 
tion is not Gaussian, although we do not expect the actual 
error to be more than a few percent. The results are given 
for the LMC and bulge (Baade's Window) in Tables 1 and 
2, respectively. 

The first rows of Tables 1 and 2 give the optical depth, 
lensing rate and average timescale, under the assumption 
that the detection efficiency e(te) = 1 for all timescales, 
whilst for the other rows detection efficiencies have been in- 
cluded in the calculation. For the LMC line of sight we use 
the MACHO efficiencies (Alcock et al. 1995), and for the 
bulge line-of-sight calculations we use the OGLE efficiency 
estimates (Udalski et al. 1994). The MACHO collaboration 
has also published efficiency estimates for its bulge obser- 
vations (Bennett et al. 1995), although they are only sam- 
pling efficiencies and do not take account of photometric 
effects. We shall therefore only compare the predictions of 
our Galactic model with the OGLE observations. 

2.3.1 Lensing towards the LMC 

The first thing to notice from Table 1 is that the combined 
optical depth of the thin disc, thick disc and bulge accounts 
for only 4% of the total optical depth for our adopted model. 
The halo is by far the most dominant source of Galactic lens- 
ing. The lensing rate is likewise dominated by the halo, with 
the combined bulge and disc contribution ranging from 1% 
for 1 Mq deflectors up to 9% for 10"* Mq lenses ( assummg 
all components are comprised of lenses of a similar mass). 
If detection efficiencies are neglected the contribution of the 
disc and bulge to the total rate is 2.4% and so including the 
efficiency estimates results in an enhancement of the rela- 
tive contribution of the bulge and disc to the overall rate for 
lenses below 0.1 Mq, and a reduction above 0.1 Mq. 

The average observed timescale is essentially given by 
the average timescale for the halo component. The aver- 
age timescale for the 3 events observed in the first year 
MACHO data is 13.8 d (Alcock et al. 1995) which, from 
Table 1 (and Section 3.1), suggests halo lenses of around 
0.01 Mq — 0.1 Mq. Taking the observed timescales as a 
guideline, it is an interesting exercise to estimate the maxi- 
mum contribution to the total rate from the bulge and disc 
components, allowing their lens masses to be different from 
those in the halo, and assuming all components are com- 
pletely comprised of compact objects. The halo rate peaks 



for m ~ 0.01 Mq, although the observed timescales are quite 
compatible with 0.1 Mq deflectors, which give only half the 
rate (around 16 events per 10'^ stars per year). The bulge 
and disc rates all peak for m ~ 10"'^ Mq lenses with a com- 
bined total of around 1.3 events per 10'^ stars per year. This 
implies a contribution of 8% to the total rate. The contribu- 
tion is as much as 13% for a dark matter-dominated bulge. 

The MACHO exposure of 9.7 x lO'' star-years (8.6 miL 
lion stars monitored over 1.1 years) actually yielded only 3 
events (Alcock et al. 1995), rather than the 16 or so events 
predicted by the above model. The relative likelihood of only 
3 events being detected, when the expectation is 16, is only 
3 X 10"*. Turning this around: given 3 events have been 
observed, the 68% confidence level (CL) lower limit on the 
expectation value for the rate is 2.9 events per year (1.4 
at the 95% CL) and the 68% CL upper limit is 4.6 (7.7 at 
the 95% CL). Fixing the combined disc and bulge rate at 
1.3 events per year, the contribution of the halo is anything 
between 55% (7%) and 72% (83%) of the total rate at the 
68% CL (95% CL). Translating this into an estimate of the 
fraction of our adopted halo comprised of compact matter, 
we find the halo fraction fu to lie between 0.10 < fu < 0.21 
(68% CL) or 0.006 < fh < 0.40 (95% CL), assuming halo 
lenses to have a mass of around 0.1 Mq. For a dark matter- 
dominated bulge the contribution of the disc and bulge is 
as much as 2.3 events per year (assuming 10"^ Mq deflec- 
tors) and so the inferred halo fraction is 0.04 < fu < 0.14 
(68% CL) or < /fe < 0.34 (95% CL). It therefore seems 
clear that the halo is not completely comprised of compact 
matter, although the 'no halo compact matter' hypothesis 
is ruled out at the 96% CL (or at the 80% CL in the case of 
a dark matter-dominated bulge). 

In fact, there a number of reasons why the combined 
disc and bulge rate would not be as high as quoted above. 
Firstly, there is evidence that sub-stellar objects (or brown 
dwarfs) probably do not dominate the density in either the 
thin disc (Kroupa, Tout & Gilmore 1993; Kerins & Carr 
1994) or the bulge (Bahcall, Flynn, Gould & Kirhakos 1994), 
implying that the bulk of these lenses must be larger than 
0.1 Mq (i.e. stellar remnants or hydrogen-burning stars) and 
thus have a lower lensing rate. In the case of hydrogen- 
burning stars, the lenses cannot be nearby, otherwise they 
would be directly observable and would not, in general, give 
rise to achromatic lensing events (since the lens and source 
would not, in general, have the same colour). Along the LMC 
line of sight, the disc and bulge density distributions peak 
close to the observer, so there is a significant contribution to 
the rate from visible lenses which in general produce non- 
achromatic light curves. For this reason we do not expect 
the visible components of the disc and bulge to give rise to 
many achromatic lensing events towards the LMC. These 
reductions may, in part, be offset by the rate due to lenses 
in the LMC itself, which is not considered here [see e.g. 
De Riijula et al. (1994)]. The LMC may make a substantial 
contribution to the rate if it has its own dark halo of com- 
pact matter, although a scenario in which the LMC halo 
is completely comprised of compact matter, but our own 
halo is not, seems contrived. In the absence of such a halo, 
the contribution of the LMC (from lenses in the LMC disc) 
would be at most comparable to our own disc. 
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Table 1. The optical depth, lensing rate and average timescale for events towards the LMC (/ = 280°, b = —33°). 
The first row shows these quantities under the assumption of perfectly efficient detection (e = 1) and a lens mass 
of 1 Mq. For other masses m these values scale as (m/MQ)"-'^/^ for F and (m/MQ)-'^/^ for (te)- The other rows 
incorporate the efficiency estimates of the MACHO collaboration (Alcock et al. 1995). r and F are displayed as 
the sum of the individual component contributions (bulge + thin disc + thick disc + halo) with their combined 
total in brackets, (tg) is similarly displayed with the rate-weighted average appearing in brackets. For the dark 
matter-dominated bulge, the bulge optical depth and lensing rate values tabulated below should be multiplied by a 
factor of 7. 



Iog(m/M0) 


t/10-8 


r/(10-^ events star"! yr~^) 




(te)/days 




0(. = 1) 


0.33 + 1.84 + 0.35 + 56 (59) 


0.08 + 0.35 + 0.09 + 21 


(21) 


96, 120, 95, 63 


(64) 





as above 


0.01 + 0.03 + 0.01 + 4.2 


(4.3) 


64, 75, 65, 46 


(46) 


-1 


as above 


0.07 + 0.29 + 0.07 + 16 


(16) 


28, 34, 28, 20 


(20) 


-2 


as above 


0.14 + 0.72 + 0.16 + 28 


(29) 


11, 13, 11, 8.0 


(8.2) 


-3 


as above 


0.17 + 0.97 + 0.18 + 24 


(25) 


4.4, 5.1, 4.2, 3.5 


(3.5) 


-4 


as above 


0.06 + 0.47 + 0.06 + 5.8 


(6.4) 


2.5, 2.6, 2.4, 2.5 


(2.5) 



Table 2. The optical depth, lensing rate and average timescale for events towards the bulge (/ = 0.°7 to l.°3, 
b = — 4.°2 to — 3.°5). The layout is analogous to Table 1 with the first row values assuming perfectly efficient 
detection and a lens mass of 1 Mq. The other rows incorporate the efficiency estimates of the OGLE collaboration 
(Udalski et al. 1994) . The values assume that all sources reside in the bulge component and have a V-band luminosity 
function dn/dC oc . For the dark matter- dominated bulge, the bulge optical depth and lensing rate values 

tabulated below should be multiplied by a factor of 7. 



Iog(m/M0) 


t/10-8 


r/(10-^ events star"! yr~^) 




(te)/days 




0(. = 1) 


5.5 + 84 + 3.5 + 15 (108) 


3.4 + 38 + 1.5 + 13 


(56) 


37, 51, 55, 27 


(45) 


1 


as above 


0.35 + 2.7 + 0.10 + 1.6 


(4.8) 


75, 92, 90, 67 


(82) 





as above 


1.4 + 16 + 0.61 + 5.0 


(23) 


36, 47, 49, 28 


(42) 


-1 


as above 


3.1 + 41 + 1.6 + 9.9 


(56) 


14, 18, 19, 10 


(16) 


-2 


as above 


4.4 + 66 + 2.6 + 12 


(85) 


5.5, 6.5, 7.3, 4.0 


(6.1) 


-3 


as above 


2.3 + 48 + 2.1 + 3.5 


(56) 


3.1, 3.1, 3.5, 2.8 


(3.1) 


-4 


as above 


0.36 + 7.0 + 0.4 + 0.39 


(8.2) 


2.7, 2.6, 2.8, 2.7 


(2.6) 



2.3.2 EROS short timescale search 

The EROS collaboration is undertaking two microlens- 
ing programs. One is a long-timescale search towards the 
LMC and SMC, which is sensitive to event timescales 
te ^ 2 days and has so far uncovered 2 LMC events 
(Aubourg et al. 1993). The other is a high time resolution 
CCD search towards the LMC, which is sensitive to events 
of 30 min ^ te ^ 7 days duration, roughly corresponding 
to lenses of mass 10"** M© < m < 10"^ M© (Aubourg 
et al. 1995). The collaboration has detected no events with 
A > 1.2 having monitored 8.2 x 10* stars for 10 months, thus 
ruling out objects with 5 x 10"** M© < m < 7 x 10"* M© 
from providing all of the halo dark matter (90% CL). Inter- 
estingly, they did find 5 candidate events with 1 < yl < 1.2, 
although these are below their adopted amplification thresh- 
old. 

Simmons, Newsam & Willis (1995) have calculated the 
flux amplification for sources with radius comparable to or 
greater than the Einstein radius of the lens (as projected 
onto the source plane). They find that, for a given am- 
plification A, there is a critical source radius = Rc — 
3 ReL/x at which the impact parameter u decreases rapidly 
as Rs increases, until the source can no longer be am- 
plified by as much as A for any value of u [see Fig. 6 
of Simmons, Newsam & Willis (1995)]. For halo lenses 



Rc ~ 30 Rq {mi/lQ-'^ MqY'^, assuming x ~ 8 kpc (which 
is a typical lens distance towards the LMC). In this case, 
sources with Rs ~ 30 Rq can be amplified by yl ~ 1.2 
for a range of impact parameters < m < 2. If we take 
LMC main-sequence stars to have a power-law mass func- 
tion dris/dms (x m~^' [where the exponent used is the 
Scalo value (Scalo 1986) for the local Solar neighbourhood, 
although determinations of LMC mass function appear to 
give similar results (Hill, Madore & Freedman 1994)], and a 
mass-radius relation oc i?^ «^ to give AusjARs oc Rg ' , 
then the fraction of source stars with a radius sufficiently 
large that an intervening lens of mass mi can only amplify 
its brightness by yl < 1.2 (irrespective of u) is 

ff"*"" dnJdR,dR, 
fiRs > Rc) = ^ . (2.9) 

For the mass function and mass-radius relation we have 
adopted, eqn (2.9) is insensitive to the maximum radius 
Rmax and so / ~ [Rcl Rmtn)~^"' , whcrc Rmtn is the radius 
corresponding to a source star of luminosity Cmm- For the 
EROS CCD survey V{lim) ~ -|-19, so ~ 50 Cq which 

implies Rmtn ~ 5 Rq. Hence, / ~ 0.3 (mi/10"'' Mq)"^/*. If 
one expects to detect N events due to lenses of mass mi 
then the probability that they all lens source stars with 
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Rs > Rc and so produce an amplification A < 1.2 is 
P{N) = ~ 0.3^ {mi/lO-" Mo)-^^/''. 

The EROS team rejects its five candidates with A < 
1.2 on the basis that, if their amplifications are given by 
the point-source approximation [eqn (2.2)], they occupy an 
improbable portion of the impact-parameter distribution. 
From our crude analysis, we see that the probability that 5 
lensing events all occur for source stars with Rs > Rc (in 
which case the point-source approximation is invalid and 
only low-amplification events are expected) is only around 
10"'^ for 10~^ Mq lenses but rapidly tends to unity for lens 
masses approaching 10"'^ Mq. This estimate is quite crude 
in many respects but does seem to indicate that at least 
some of the 5 EROS candidates could be interpreted as 
true microlensing events if the lenses have a low mass, say 
mi ~ 10"'' Mq, and are lensing source stars with large ra- 
dius (Rs ~ 10 -Rq). The fact that the source stars all appear 
to be much brighter than average is also consistent with this 
idea. 

In fact, only 30% of the stars monitored by EROS are 
main sequence, the rest are giants. Our estimate may there- 
fore be on the pessimistic side. Intriguingly, the number of 
events due to 10~ Mq lenses expected by EROS is 5.6 if 
the halo is completely comprised of such objects! 



2.3.3 Lensing towards the Galactic bulge 

Table 2 shows the expected optical depth, lensing rates and 
average event times for the bulge line of sight. The actual re- 
gion of the bulge corresponds to the Baade's Window area 
monitored by the OGLE collaboration, in which 7 of the 
10 events reported by Udalski et al. (1994) were found. In 
the case of a non dark matter-dominated bulge, the thin 
disc provides the dominant contribution to the optical depth 
(78%), although the halo is also an important source (con- 
tributing 14%). The disc likewise dominates the rate, pro- 
viding up to 85% of the total rate in the case of 10~* Mq 
lenses. The disc dominates the rate even for the case of a 
dark matter-dominated bulge, although the bulge contribu- 
tion is comparable. The effect of the efficiencies is to en- 
hance the rate of the halo and bulge, relative to the thin 
disc, for masses above 1 Mq but reduce their contribution 
below 1 Mq. The opposite is true for the thick disc, the 
contribution of which increases relative to the thin disc for 
lower mass lenses. However, since the thick disc is not a ma- 
jor contributor to the overall rate, the thin disc provides a 
greater proportion of the observed total rate for lower mass 
lenses. 

The average event timescale for the 7 OGLE lenses 
found in the Baade's Window region is 33 d, implying typ- 
ical lens masses in the 0.1 Mq — 1 Mq range from Table 2 
(see also Section 3.1). Of these, 2 have minimum impact pa- 
rameters larger than Re (u > 1), leaving 5 events with an 
average timescale of 31 d. On the basis of the average ob- 
served timescale, the expected rate is ~ (2.3 — 5.6) x 10~^ 
events per star per year. The OGLE collaboration monitored 
9.5 X 10^ stars in 1992 (which yielded 3 of the 5 events) and 
7.5 X 10^ in 1993 (allowing for field overlaps and contami- 
nation due to disc sources) and so the expectation rate is 
2.2 - 5.3 events in 1992 and 1.7 - 4.2 events in 1993. These 
rate intervals include the observed values and so our stan- 



dard Galaxy is not actually at odds with the OGLE data as 
far as the observed rate is concerned. 

However, there does appear to be a discrepancy be- 
tween the model and the observations concerning the op- 
tical depth. The total optical depth for the model is only 
(1.1 — 1.4) X 10~^, where as the value of r inferred by 
the OGLE collaboration (based on all of its events) is 
(3.3 ± 1.2) X lO"'' (Udalski et al. 1994) [a similar value has 
been obtained by the MACHO collaboration, based on its 
sample of 13 'clump giant' lensing events (Bennett et al. 
1995)]. It is not clear how serious this discrepancy is, since 
one might not be surprised by an optical depth estimate, 
inferred from only 10 events, deviating significantly from 
the expectation value. It is unclear how large a deviation 
is allowed (since it depends on the relative contributions of 
the various components), although it should be much larger 
than suggested by Poisson estimates. 

It is interesting to note from Table 2 that, whilst the 
(non dark matter-dominated) bulge contributes less than 
6% to the total rate (for a lens mass m < 1 Mq), the rate- 
weighted average timescale for all components is typically 
fairly close to the average timescale for the bulge. Therefore, 
in the case of a dark matter-dominated bulge 0.1 Mq — 
1 Mq lenses are still favoured. However, for this mass range, 
bulge lenses become a major contributor to the rate with 
(1.0-2.2) X lO"'' events per star per year, and the predicted 
total rate for the OGLE 1992 and 1993 observing seasons 
becomes 3.0 — 7.0 and 2.3 — 5.6 events, respectively. Note, 
however, that recent HST observations of the local bulge 
population appear to reduce the possibility of much bulge 
dark matter (Bahcall, Flynn, Gould & Kirhakos 1994). We 
have seen that the observed rate towards the LMC implies 
that the halo fraction is likely to be ~ 0.2 rather than 
1. This would result in a small reduction in the overall rate 
for lenses between 0.1 — 1 Mq but the model would still be 
quite consistent with the observed rate towards the bulge. 
Of course, one would also expect a small contribution to the 
rate from the lensing of disc sources. 

Alternative Galactic models invoked to reconcile the 
predicted and observed optical depth include a bar-like bulge 
(Paczyhski et al. 1994a; Evans 1994; Zhao, Spergel & Rich 
1995), for which there is some independent observational 
support (Blitz & Spergel 1991; Weiland et al. 1994). A more 
radical suggestion for increasing the optical depth involves 
assuming a maximal disc (Gould 1994a; Alcock et al. 1994). 
In this scenario, almost all of the rotation curve velocity in- 
side Ro = 8 kpc is provided by the disc, thereby allowing a 
much larger density (and implying a substantially reduced 
halo). However, Paczyhski et al. (1994b) claim that there is 
evidence of a dramatic decrease in the disc density within 
the inner ~ 4 kpc of the Galaxy. Such a hollow disc could not 
provide the required optical depth and so a bar-like bulge 
aligned close to the observer's line of sight would appear the 
more plausible alternative. (A hollow disc would likewise 
lower the rate and optical depth estimates calculated in this 
paper, as well as altering the average timescale estimates.) 
These models ought to be testable in the near future with 
increased angular coverage by the microlensing experiments. 
Whilst there is other evidence which points towards a Galac- 
tic model which differs from the standard one investigated 
here, we believe that, on the basis of microlensing observa- 
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Figure 1. The number of simulated halo events of mass m with 
event duration tg along the (a) LMC and (b) Galactic bulge 
(Baade's Window) lines of sight. The adopted values for tg re- 
flect the event durations that are typically being observed. The 
most likely lens mass is 0.05 Mq [te/15 d)^ for the LMC direction 
and 1.2 Mq (te/30 d)^ for the bulge. 
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Figure 2. The number of simulated thin-disc events of mass m 
with event duration tg along the (a) LMC and (b) Galactic bulge 
(Baade's Window) lines of sight. The most likely lens mass is 
0.02 Mq (te/15 d)2 for the LMC direction and 0.3 Mq (te/30 d)^ 
for the bulge. 



tions alone, the case for a radical change in the model is not 
yet an overwhelming one. 



3 LENS MASS ESTIMATES AND EVENT 
AMPLIFICATIONS 

Having looked at the optical depth, lensing rate and average 
event duration, we now turn our attention to the individual 
lenses themselves. 



3.1 Lens mass 

Mass estimates of individual lenses cannot be made very 
precisely without additional parallax (Gould 1994b), source- 
polarisation (Simmons, Newsam & Willis 1995) or a combi- 
nation of high-precision astrometric and photometric obser- 
vations (H0g, Novikov & Polnarev 1995). This is because 
the event duration not only depends on the lens mass but 
also on its relative tangential velocity, as well as the lens and 
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Figure 3. The number of simulated thicli-disc events of mass m 
with event duration tg along the (a) LMC and (b) Galactic bulge 
(Baade's Window) lines of sight. The most likely lens mass is 
0.03 Mq (te/15 d)2 for the LMC direction and 0.3 Mq (te/30 d)^ 
for the bulge. 



Figure 4. The number of simulated bulge events of mass m with 
event duration tg along the (a) LMC and (b) Galactic bulge 
(Baade's Window) lines of sight. The most likely lens mass is 
0.03 Mq {te/15 d)2 for the LMC direction and 0.5 Mq (te/30 d)^ 
for the bulge. 



source distances. The lens mass is given by 

^ = 0.08 Me f^V ( " ( S , 
"V20d; 1^10 kpcy V220 kms-V 

(3.1) 

where D = x{L — x)/L. The uncertainty is contained in the 
unknown lens distance x and tangential motion Vt (and also, 
for observations towards the bulge, the source distance L). 

Griest (1991) has shown how an estimate of the lens 
mass can be made from the variation of dT/dte with mass 



for a given value of te. De Riijula, Jetzer & Masso (1991) 
have shown how the moments of the event duration can 
also be used to extract information on the lens mass func- 
tion. This method has been used by Jetzer (1994a; 1994b) 
to estimate the mass of the first LMC and bulge events. In 
this section, we wish to calculate the lens mass probabil- 
ity distributions for each of the four Galactic components, 
both towards the LMC and Galactic bulge. In the previous 
section we used the average observed timescales to indicate 
likely lens mass ranges and to compare observed and pre- 
dicted lensing rates. Here, we shall see how well those very 
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approximate mass range estimates correlate with more rig- 
orous statistical determinations. There are numerous theo- 
retical and observational constraints on many of the bary- 
onic dark matter candidates (Carr 1994), and these can be 
used in conjunction with the lens mass estimates to try to 
determine the dominant lensing components. 

For our lens mass estimates we employ a Monte-Carlo 
simulation of 2 x 10^ lensing events for each component and 
each line of sight. In the case of the bulge line of sight a 
source distance L is selected at random for each event, as- 
suming the distribution of eqn (2.8), with /? = 2 and ps 
given by eqn (2.7). For the LMC line of sight L is taken to 
be fixed at 50 kpc. For the lens distance x a random value 
between and L is selected assuming a probability distri- 
bution P{x)dx = dT{x)/T{< L) oc D{x, L)pi{x), where pi 
depends upon the component being considered. A random 
velocity is then assigned to the lens (and also the source in 
the case of the bulge line of sight) using eqn (2.5) and the 
contribution of any systematic motion (i.e. any rotational 
velocity component of the lens or motion of the line of sight 
at the lens position) taken into account to give a relative 
transverse lens velocity Vt. For a given timescale the lens 
mass is then specified by eqn (3.1). For the LMC line of sight 
we adopt te = 15 d and for the bulge direction = 30 d as 
being representative of the observed event durations. 

The results are shown in Figs. (1-4) for halo, thin disc, 
thick disc and bulge lenses, respectively. For the LMC direc- 
tion [Figs. l(a)-4(a)] the most likely lens mass is 0.05 Mq, 
0.02 Mq, 0.03 Mq and 0.03 M© for the halo, thin disc, thick 
disc and bulge components, respectively. The implication is 
that Galactic lenses observed towards the LMC are likely 
to be brown dwarfs (no matter which component is respon- 
sible for the observed events), although the 95% CL mass 
range is around 2.5 orders of magnitude for the halo and 2 
orders of magnitude for the other components, so microlens- 
ing observations alone do not rule out M-dwarf stars. How- 
ever, high-latitude searches for M-dwarfs in the optical and 
near-infrared (Bahcall, Flynn, Gould & Kirhakos 1994; Hu, 
Huang, Gilmore & Cowie 1994) have yielded significantly 
fewer stars than required to explain the observed optical 
depth. It therefore seems that the LMC microlensing events 
are due to brown dwarfs. 

In the case of the lensing events observed towards 
Baade's Window [Figs. l(b)-4(b)], the most likely lens mass 
is 1.2 Mq, 0.3 Mq, 0.3 M© and 0.5 M© for the halo, thin 
disc, thick disc and bulge components, with a 95% CL range 
of around 2.5 orders of magnitude in each case. In the case 
of the halo, a lens mass of 1.2 Mq implies either hydrogen- 
burning stars or stellar remnants. The idea of the dark halo 
being comprised of ordinary main-sequence stars is clearly 
at odds with observations, whilst a halo comprised of stel- 
lar remnants would over-enrich the interstellar medium with 
metals, unless the mass function is finely tuned to around 
5 Mq (Ryu, Olive & Silk 1990). In any case, the disagree- 
ment between the halo lens masses as inferred along the 
LMC and bulge lines of sight implies that the halo is not 
dominating the rate along at least one of these directions. 
Timescales of around 10 d towards the bulge are required 
to provide agreement between the inferred halo lens masses 
along both lines of sight, and to provide consistency with 
other observational and theoretical constraints. Since the 
observed timescales are typically higher than this, the halo 



cannot be dominating the rate towards the bulge. The ob- 
served rate is anyway much higher than can be accommo- 
dated by the halo. For the other components M-dwarf or 
brown dwarf mass scales are preferred. However, optical and 
near-infrared searches for M-dwarfs imply that they cannot 
be major contributors to the lensing rate (Bahcall, Flynn, 
Gould & Kirhakos 1994; Hu, Huang, Gilmore & Cowie 
1994), unless the local mass function is unrepresentative of 
the mass function averaged along the line of sight towards 
Baade's Window. Similarly, the disc mass function must 
rise sharply below the hydrogen-burning limit (0.08 Mq) 
if brown dwarfs are to provide the required lensing rate in 
our model (Kroupa, Tout & Gilmore 1993; Kerins & Carr 
1994). 



3.2 Event amplifications 

For a peak amplification A the minimum line-of-sight dis- 
tance of the source star from the lens is smaller than Re 
by a factor n = 2^'^[A(A^ - l)'^'^ - if'^ [the inverse of 
eqn (2.2)]. Hence, one can infer the impact parameter di- 
rectly from observations. 

The probability of a lensing event having an impact 
parameter in the range (u, u + du) is P{u) du oc du and so for 
a large number of events one would expect the mean impact 
parameter to be (m) = 0.5 or A({u)) = 2.18 for a maximum 
'threshold' impact parameter Ut = 1. The probability of 
observing n independent events with impact parameters 
and mean ¥ = Ui/n relative to n events with a mean 

(m) = Ut /2 is 



P-/P(u) 



(«/(«))"-' 
[{ut-u}/(u)r 



(0 < ¥< (m)) 
((m) < ¥ < Mt) 



(3.2) 



Thus, an improbable sample would have P—/P(u) ^ 1- The 
above likelihood estimator holds true provided there are no 
observational selection effects which favour some amplifica- 
tions more than others. Inevitably, however, such selection 
effects do arise and so eqn (3.2) is as much an estimate of 
the importance of selection effects as it is a measure of the 
likelihood of the sample. 

The MACHO collaboration adopts a minimum amplifi- 
cation cutoff oi A = 1.5, implying Ut = 0.83 and (u) = 0.42. 
The 27 bulge events discussed by Sutherland et al. (1995) 
have a mean impact parameter u ~ 0.39 [from Fig. 4 of 
Sutherland et al. (1995)]. The relative likelihood of this sam- 
ple is therefore P—/P(u) — 0-3. This likelihood is extremely 
sensitive to u for a sample as large as 27 events and so 
the fact that P—/P(u) is anywhere close to 1 is an indi- 
cation of a very uniform sample. The MACHO team has 
detected a further 18 candidates (Bennett et al. 1995) but 
the amplifications for these await publication. MACHO has 
also detected 4 LMC events with amplifications 1.52, 1.99, 
2.98 and 7.2, giving a mean impact parameter u = 0.47 and 
P—/P(u) = 0.67, which again implies a pretty uniform sam- 
ple. The combined likelihood for the 4 LMC and 27 bulge 
events is P—/P(u) — OA. 

The EROS team adopts Ut = 1.0, so for its sample 
(u) = 0.5. It has so far detected only 2 lensing events from 
their longer timescale search and these have amplifications 
of 2.7 and 3.0 giving u = 0.37 and P—/P(u) = 0.75. This high 
likelihood, more than anything, reflects the fact that we are 
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dealing with only 2 events. EROS finds no candidate events 
with A > 1.2 from its short timescale CCD observations 
(Aubourg et al. 1995) although, as discussed in Section 2.3.2, 
the 5 light curves with 1 < yl < 1.2 could be interpreted as 
very low mass objects lensing large LMC stars. In this case, 
the impact parameter distribution would not obey eqn (3.2) 
(Simmons, Newsam & Willis 1995). 

The OGLE team has not adopted an amplification 
threshold and so very low amplification events are accepted 
provided their light curves have sufficiently small photomet- 
ric errors. However, in the absence of selection effects, the 
impact parameter distribution should still be uniform below 
some adopted threshold Ut. There are now 12 OGLE bulge 
events, one of which is thought to be a double lens. Of the 
remaining 11, 3 have A < 1.34 and so taking an adopted 
threshold of Mt = 1 leaves 8 events. The mean impact pa- 
rameter for these 8 is ¥ = 0.35, implying P—/P(u) = 0-08 
from eqn (3.2). This appears to imply that the OGLE sam- 
ple is quite improbable. However, OGLE has conducted a 
more sophisticated analysis of its impact parameter distri- 
bution (Udalski et al. 1994) which takes account of photo- 
metric selection effects. OGLE determines Ut for each event, 
based on the photometric errors at the time of observation, 
and takes u/ut, rather than u, as an indicator of the uni- 
formity of the sample, u/ut should be uniformly distributed 
between and 1 and so (u/ut) should converge to 0.5 for a 
sufficiently large sample. Eqn (3.2) can then be used to test 
the distribution of u/ut rather than u. For the 9 single-lens 
events reported by Udalski et al. (1994), u/ut = 0.53 and 
so P^i^^ l P(ulut) = 0.58, which indicates that the sample is 
indeed consistent with the microlensing hypothesis. 



4 CONCLUSIONS 

We have investigated the expected lensing properties of a 
'standard Galaxy' model with spherically-symmetric halo 
and bulge/spheroid components. The simplifying assump- 
tion of spherical-symmetry is naive and, in the case of the 
bulge, is not favoured by direct observations, but we are here 
interested in the extent to which microlensing rules out such 
a model. 

Using the observed event durations towards the LMC 
we estimate the total rate due to dark matter to be around 
16 events per year for MACHO observations, most of which 
comes from the dark halo. MACHO detected only 3 events 
in its first year and so, if one takes the contribution of the 
disc and bulge components to be fixed, this limits the halo 
fraction in compact objects to fh < 0.4 at the 95% CL. This 
limit is fh < 0.34 if the bulge is dominated by dark matter. 
However, a 'no halo compact matter' hypothesis is ruled 
out at the 80% CL or greater unless one assumes that the 
LMC itself has a substantial halo comprised of such objects 
(De Riijula, Giudice, MoUerach & Roulet 1994). A scenario 
in which the LMC halo is comprised of compact objects, but 
our own Galactic halo is not, seems rather contrived. 

We have also computed lens-mass probability distribu- 
tions for each component, based on the observed event du- 
rations, and have found that, for the LMC direction, the 
lens is likely to be a brown dwarf, irrespective of which 
component the lens actually resides. Whilst M-dwarfs are 
not ruled out by microlensing observations, their contribu- 



tion is strongly constrained by near-infrared searches (Bah- 
call, Flynn, Gould & Kirhakos 1994; Hu, Huang, Gilmore & 
Cowie 1994). 

For the bulge line of sight we find an expected lensing 
rate of around 2 — 5 events per year for OGLE observations 
(again using the observed timescales as a guide and count- 
ing as events only those detections for which m < 1), which 
is comparable to what they are observing. The predicted 
optical depth is a factor 2 or 3 less than the current ob- 
servational determinations but these estimates are based on 
only 10 events (or 13 in the case of MACHO) and so this 
discrepancy is probably not too serious. More detections are 
required to conclusively show up any discrepancy between 
spherically-symmetric and asymmetric bulge models. 

Lensing mass determinations for the bulge line of sight 
point to M-dwarf or brown dwarf mass scales in the case 
of the disc or bulge. To provide consistency with near- 
infrared limits one needs to assume that the lenses are brown 
dwarfs rather than M-dwarfs and, even then, one requires a 
sharp turn-up in the slope of the mass function below the 
hydrogen-burning limit to provide the required rate, unless 
the local mass function is unrepresentative of the mass func- 
tion averaged along the bulge line of sight. For lenses in the 
halo component a mass scale around 1 Mq is preferred on 
the basis of the observed event durations. Such objects could 
not be hydrogen-burning and so must necessarily be stel- 
lar remnants. However, the possible contribution of stellar 
remnants is already strongly constrained by ISM metallicity 
observations and so the implication is that the halo cannot 
be dominating the rate towards the bulge. This is consistent 
with a halo only partially comprised of compact matter, as 
has already been inferred from LMC microlensing observa- 
tions. 

A naive analysis of the distribution of impact parame- 
ters for the EROS, MACHO and OGLE candidates, which 
ignores possible photometric bias, shows that the EROS and 
MACHO candidates are quite consistent with the microlens- 
ing hypothesis. The OGLE sample is also consistent with 
microlensing provided careful account is taken of photomet- 
ric effects. This is primarily because the OGLE collaboration 
has a different selection procedure from MACHO and EROS 
in that it does not adopt an amplification cutoff. 

The EROS team has recently placed strong limits on 
the density contribution of very low mass (10"** Mq < m < 
10"^ Mq) objects on the basis that its short timescale CCD 
search has uncovered no lensing events with A > 1.2. How- 
ever, EROS does have 5 low-amplification light curves which 
might be best explained as the microlensing of large source 
stars by very low mass lenses. In such a case the point- 
source approximation breaks down and so one would expect 
only low-amplification events to occur. On the basis of re- 
cent work on finite-size source lensing by Simmons, Newsam 
& Willis (1995) we estimate that such low-amplification 
events are extremely common for lens masses approaching 
10"'' Mq. Follow-up observations of the source stars them- 
selves should help to discriminate between this suggestion 
and the possibility of intrinsic variability. 
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